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ABSTRACT 

We discuss a model of an X-ray illuminating source above an accretion disk of a rotating black hole. 
Within the so called lamp-post scheme we compute the expected (observed) polarization properties 
of the radiation reaching an observer. We explore the dependencies on model parameters, employing 
Monte Carlo radiation transfer computations of the X-ray reflection on the accretion disk and taking 
general relativity effects into account. In particular, we discuss the role of the black hole spin, of the 
observer viewing angle, and of the primary X-ray source distance from the black hole. We give several 
examples of the resulting polarization degree for two types of exemplary objects - active galactic 
nuclei and Galactic black holes. In order to assess potential observability of the polarization features, 
we assume the sensitivity of the proposed New Hard X-ray Mission (NHXM). 

We examine the energy range from several keV to ~ 50keV, so the iron-line complex and the 
Compton hump are included in our model spectra. We find the resultant polarization degree to 
increase at the higher end of the studied energy band, i.e. at > 20 keV. Thus, the best results for 
polarimetry of reflection spectra should be achieved at the Compton hump energy region. We also 
obtain higher polarization degree for large spin values of the black hole, small heights of the primary 
source, and low inclination angles of the observer. 

Subject headings: polarization — relativistic processes — X-rays: galaxies — X-rays: binaries 



1. INTRODUCTION 

In active galactic nuclei (AGN) and Galactic black 
holes (G BH), X-rays are p roduced near the black hole 
horizon (jFrank et al.l 12002 ). The spectrum of these ob- 
jects typically consists of several constituents, namely, 
a multi-color thermal component arising from an opti- 
cally thick accretion disk and a power-law component, 
thought to origi nate in an optically th i n medium, a .k.a. 
a 'corona' (e.g., IBlandford etHlfToOOt IKrolikI [TOOl . A 
plausible scheme suggests that the latter component is 
produced in a hot Comptonizing layer above the disk. 
Various spectral features originate by the mutual inter- 
action of the accretion disk medium and the corona. In 
particular, very prominent are the reflection features of 
iron, which occur in the 6~7 keV range in the local frame 
of the emitti ng material (E,emillard & McClintock 200l; 
iFabianI [20081 ). 

The mutual proportions among these components in 
a particular object depend on the physical (spectral) 
state of the system, as well as the intervening medium 
along the line of sight. These proportions are also influ- 
enced by strong gravity near the black hole, where light 
bending effects and energy shifts can play a significant 
role. Disentangling the observed spectrum into its in- 
dividual components is not a straightforward task, as it 
poses a degenerate problem and does not generally allow 
us to constrain all parameters. Nevertheless, accurate 
determination of the spectral components is an essen- 
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tial step in measuring the black hole spin, which is one 
of the most important challenges in the present-day as - 
trophysics of accreting black holes (jFabian et al.|[2004[l . 
It has been proposed that X-ray polarimetry could pro- 
vide additional constraints on the accretion disk models, 
limit the range of acceptable parameter values, and help 
us reducing the ambiguities that cannot be resolved by 
spectroscopv alone (Dovciak et al.l l2008t iLi et al.l 120091 : 
iSchnittman fc Krolikll2009i r 

In the inner regions of black hole accretion disks, 
eme rging photons are strongly influenced by strong grav- 
ity (jlCato et al.l 119981 ). The framework of General Rel- 
ativity (GR) is n eeded to properly determine the re- 
sulting spectruni (jFabian et aLlhoOOt iKarasI 120061: iMatti 
l2006HMilleij|2007D . Relativistic effects include the energy 
shifts, both due to special relativity (Doppler effect) and 
GR (gravitational redshift), as well as the light bending 
aberration effects that are particularly prominent near 
the photon orbit: Tph = SGM/c^ = 4.43M/MQkm for 
a non-rotating black hole (the radius of rph decreases 
with the black hole spin increasing). The relevance of 
GR effects is even more pertinent for the calculations of 
the observed polarimetric quantities (Connors & StarH 
Tpr?; 'Connors et al."'1980; Laor et al. 1990; Matt 1993 
Ago! fc Krolik 2000; Dovciak et al. 2004b, 2008; Li et aT| 
20091 : ISchnittman fc Krolikll2009l 120101 . This is mainly 



due to sensitivity of the polarimetry parameters to the 
geo metrical arrangement of the source. 

In iDovciak et all (|2004bl l. we studied the GR effects 
on the polarization properties of the radiation reflected 
by an accretion disk near a rotating (Kerr) black hole. 
We considered illumination of the disk surface by a hot 
corona. The corona acts, via the inverse Compton effect, 
on thermal disk photons and it produces as an outcome 
the above-mentioned power-law component. 

In this context, a special scheme has been dubbed the 
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'lamp - post model' (iMatt et a^lQQlHMartocchia fc Matti 
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Miniutti fc Fabiaijl20o3 iDovciak et al.ll2004bD . Accord 
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ing to this model, the illumination is caused by a point- 
like source located on the rotation axis of the black 
hole, which is also the accretion disk axis. This kind 
of a source can be identified, e.g., with the base of 
a jet, although the exact interpretation of the 'lamp' 
does not need to be specified at this stage of the phe- 
nomenological scenario. Generalizations of the lamp- 
post model to the case of extended sources, off-axis il- 
lumination, and non-planar disks have been also investi- 
gated (jPabrowski fc Lasenbyll200"TI : iNiedzwiecki fc Zvckil 
l2008f ). Different, but also highly relevant and interest- 
ing geometry of an extended corona has been recently 
studied by Schnittman & Krolik (2010), taking into ac- 
count relativistic effects and even the self-irradiation of 
an accretion disk by its own returning radiation. 

The motivation for the present paper arises from the 
advent of X-ray po larimetry in the field of AGN and 
GBH (Bel lazzini et al. 2006; BcU azzini fc Span drc 2010). 
The paper revisits the light bending scenario of accret- 
ing black holes and it improves the calculations of the 
flux and polarization of the reflected disk emission, with 
the emphasis on the observability of GR effects and the 
possibility to constrain the model parameters with fu- 
ture X-ray polarimeters on-board satellites, such as the 
currently proposed New Har d X-ray imaging and po lari- 
metric Mission (NHXM) (Tag haferri et al.l l2010a llbl). or 
the confirmed Gr avity and Extreme Magnetism SMEX 
(GEMS) mission (|Swank et al.l[2010l) . 

The paper is organized as follows. In section [2] we 
briefly introduce the lamp-post scheme for black-hole 
accretion disks. In section |3] we describe the primary 
component which acts as a source of the X-ray irradi- 
ation. Here we also provide the relevant equations and 
we calculate the illumination of the disk. In section |4] 
we describe details of the radiation transfer computation 
that we apply in order to derive the reflection compo- 
nent. These considerations are then employed in section 
[SI where we discuss a consistent procedure connecting 
the primary source position on the rotation axis near the 
horizon with the signal observed at arbitrary direction 
far from the black hole. Results for the predicted polar- 
ization properties are then described in section |6l taking 
into account the continuum together with the iron line 
features. Finally, in section [7] we describe the observa- 
tional prospects with the currently planned X-ray po- 
larimetric missions. We demonstrate the observability of 
GR polarimctric features by modeling the expected data 
from two proto-typical sources and assuming the sensi- 
tivity of the NHXM polarimeter as an example. Section 
[5] summarizes the main results and concludes the paper. 

2. THE LAMP-POST MODEL FOR BLACK-HOLE 
ACCRETION DISKS 

The larup-post mo del was originally pro posed by 
IMatt et all (|1991l ) and lPetrucci fc Henril (|1997f ) with the 
aim to provide a simple common scheme for the origin 
of the X-ray power-law continuum and the relativistic 
spectral features seen in accreting black hole sources. A 
natural motivation for this scenario arises from the pre- 
sumed co-existence of gaseous medium in different states 
in the inner regions of the accretion disks, and hence 



the action of different radiation processes that shape the 
resulting spectra. This relates, in particular, to the spec- 
trum expected from a geometrically thin, optically thick 
accretion disk, surrounded by a hot, diluted corona, pos- 
sibly with an outgoing jet in some objects. 

The light-bend i ng s cenario was applied by 
iMiniutti fc FabianI ()2004[ ) to explain the time be- 
havior of the primary and the reprocessed emission of 
the proto-typical Seyfert 1 galaxy, MCG-6-30-15, where 
the spectral features from the reprocessing vary less than 
the primary X-ray continuum. These authors notice 
that, if the source of primary illumination is located very 
close to the black hole, the observed spectral properties 
may result as consequence of rearrangements of the 
geometrical proportions of the object. This can be 
interpreted in terms of a point-like source, representing, 
for example, shocks in the base of a jet, or the location 
where int ense flares arise via magnetic reconnection 
processes. iRossi et all ()2005D showed, using RXTE data 
for the X-ray transient XTE J1650-500, that the light 
bending scheme may well apply also to Galactic black 
hole systems. 

observer 



corona 





accretion disk 



Figure 1. Scheme of the modeL An elevated primary source 
situated on the axis of the black hole at height h illuminates the 
corotating accretion disk. The observer at infinity receives both 
the unpolarized power-law radiation from the lamp-post as well as 
the polarized reflected light from the disk. The incident angle, (5;, 
between the incoming light ray and the disk normal, the emission 
angle, 5c, between the outgoing light ray and the disk normal, 
and the relative azimuthal angle, A(f>, between the incident and 
emitted light rays define the scattering geometry in the local frame 
co-moving with the accretion disk. The normal to the disk and 
incident and emitted light rays projected onto the disk surface are 
depicted by the dashed lines. 



According to the original, most simple version of this 
model, the primary 'lamp' moves along the black hole 
axis, typically at heights of only ^ 10-20 gravitational 
radii above the disk plane. The directly seen contin- 
uum component then changes its flux very prominently, 
mainly due to the effect of changing the gravitational 
redshift. On the other hand, the reprocessing spectral 
features are less variable, therefore reproducing the ob- 
served pattern, which otherwise is difficult to apprehend. 
This scheme is attractive because it allows us to connect 
the black hole parameters to the observed spectral varia- 
tions, and so it has been pursued in a number of papers. 
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The above-mentioned interpretation, however, imposes 
the right proportions between the source height and the 
accretion disk size, so that the hght-bending effect op- 
erates efficiently in the desired way. This may not be 
generally satisfied, or it may be unlikely to happen in 
majority of sources. The check can be performed by 
comparing results of various approaches, such as tim- 
ing, broad-band spectra, and polarimetry. Recently, 
the rms variability of MCG-6-3 0-15 was investigated by 
iNiedzwiecki fc Mivakawal ()2010[ ) and the Suzaku broad- 
band spectra were discussed by iZvcki et all (|2010( ). 
These authors indeed find various complications that ap- 
pear when the simple lamp-post scenario is applied to 
complex data. 

The basic definitio n of the model set-up follows 
IDovciak et all ()2004b( ). see Figure [TJ A point-like emit- 
ting source on the black hole axis is located at the height 
h. The primary emission that is assumed to be unpo- 
larized illuminates the accretion disk where it is repro- 
cessed and part of it is emitted towards the observer. In 
the lamp post scenario the path of the light that trav- 
els from the lamp-post through some point in the disk 
up to the observer is set for each position on the disk 
unambiguously, thus defining the scattering geometry at 
every point in the disk plane, see Figure [TJ In this paper 
we consider only direct photons, i.e., we assume optically 
thick disk and we neglect the photons that might bend 
near the photon orbit in such a way that they would cir- 
cle around the black hole and only then either reach the 
observer or strike the surface of the disk. 

We introduce a number of improvement s in addition to 
our pr evious polarization computations (jPovciak et al.l 
l2004bD : 

• The reflection emission is computed with a more so- 
phisticated multi-scattering Monte Carlo code for 
the radiation transfer in the accretion disk. As a 
result of this modeling, 

— the azimuthal dependence of the emission is 
taken into account (i.e., the reflection depen- 
dence on the azimuthal directions between in- 
cident and emitted light rays); 

— the reflection is computed also for higher ener- 
gies (up to 100 keV locally measured energy) 
and the resulting Compton hump is included 
in these new computations. 

• We have further refined the computation of light 
geodesies from the axis to the disk plane, and fur- 
ther to the observer, so additional improvements 
as opposed to our previous results were achieved, 
especially 

— self-consistent computations for different spin 
values (a = 0.9987 in our previous work); 

— an extended range of possible heights of the 
primary source on the axis (we are able to 
trace the dependence on the height more pre- 
cisely); 

— geodesies from the axis to the observer are 
treated in a proper way, i.e., the lensing effect 
is included (however, we find that the previous 



approximation was already very good and this 
point does not lead to much different results). 

• We also include a contribution to the total signal 
from unpolarized fluorescent Ka and K/3 lines. 

We calculated the flux and polarization of the repro- 
cessed emission, as well as the flux of the direct emission, 
in four energy bands (2-6 keV, 6-10 keV, 10-20 keV, 20- 
50 keV) for three different inclination angles (30, 60, and 
80 degrees, respectively), and for various values of the 
height h. Two different values of the angular momentum 
of the black hole have been explored: a static, a=0, and 
an extremely rotating black hole, a=l, where, as custom- 
ary, a indicates the angular momentum per unit mass in 
units of the gravitational radius. We use these units for 
the remainder of this paper. The accretion disk in our 
computations extends from the marginally stable orbit 
(nn = ?'ms) up to the radius rout — 1000 GM/c^. 

3. PRIMARY SOURCE, OBSERVED PRIMARY RADIATION, 
AND ILLUMINATION OF ACCRETION DISK 

We assume the primary source of unpolarized radia- 
tion to be a point-like patch of the corona located at 
height h < 100 above the black hole (measured from the 
center). It radiates isotropically with the specific inten- 
sity being the usual power-law dependence on the energy, 
Iiamp{E) = Niamp E~" . In our computations we use the 
value of the power-law index a = 1 and the normalization 
factor A^iamp to be unity in the local static frame. 

The primary emission received by the observer gets 
diluted by the relativistic effects that act on photons on 
their path, mainly near the central compact body. The 
observed intensity, /prim(^^) = dE / {dt dil dv) , will be 
again a power-law in the observed energy 

Iprini{E) — A^prim E^"' , (1) 

where 

AT — "-^^'amp _ 

JVprim - 5prim " 

Q + l 

_ /, _ 2/l \ ^ sin 6>iamp d6>iamp 

^\ h^ + ay sin^o d0o ' ^ ' 

with ^lamp and 0o being the local emission angle and the 
inclination angle of the observer, respectively, both mea- 
sured with respect to the system axis. We denoted the 
energy shift (from E'lamp to i?piim) of the primary pho- 
tons detected by the observer by gprim = ^'primZ-E'iamp = 

1 /2 

^/UL^pih,a) = -gtl {h,a). Here, J/j^^p is the four- 
velocity of the motionless lamp-post and g^^ is the space- 
time metric. In Figure [2] we show the dependence of 
the factor Aprim and intensity /prim integrated in sev- 
eral energy bands on the height of the lamp-post for 
Schwarzschild and extreme Kerr black holes. One can 
see that the relativistic effects for direct radiation be- 
come small for heights ft- > 100 and the difference be- 
tween extremely rotating and non-rotating black holes 
are well visible only for heights h<Q. 

Due to the lamp-post geometry and relativistic effects, 
the primary source illuminates the accretion disk un- 
evenly. The reflected flux from the disk is proportional to 
the incident flux on the disk. The incident flux is also the 
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Figure 2. Primary emission. Left: The emission received by the observer from the primary source is diminished by the factor A'^prim 
because of the strong gravity of the central black hole. The dependence on the height of the illuminating corona, h, is shown, if the intensity 
of the primary power-law emission is E~^. This factor is almost independent on the inclination of the observer (the curves for inclinations 
1° and 85° would be separated by less than the width of the shown curves). Right: The primary component of the intensity far from the 
source (at radial infinity), integrated in the energy ranges 2-6, 6-10, 10—20, and 20-50 keV (curves with higher to lower normalization). 
The dependence on the height of the primary source is shown with the solid curves for the non-rotating Schwarzschild black hole (a=0) 
and with the dotted curves for the maximally rotating Kerr black hole (a=l). 



power-law with the same index a and with the fohowing 
normahzation factor (see lDovciak et aT]|2004blFI 

AT — ^'"c dr^iamp _ 
^^inc ~ t - 

Sin^lamp d^iamp .o\ 

X ^ —-^ . 3 

r dr 

Here, the energy shift ginc (from £^iamp to i?inc) of the pri- 
mary photons incident on the accretion disk (measured 
in local frames) is ginc = -Binc/-Biamp = C^lisk/^ilmp ^i^h 
^disk being the Keplerian four-velocity of the disk. The 
radial distance r is the Boyer-Lindquist radial coordinate 
(we use it here because eventually we integrate in Boyer- 
Lindquist coordinates on the disk and not in the local 
frame of the disk). We show the ratio of the incident 
radiation to that received by the observer directly from 
the lamp-post in Figure [3l 

Compared to the Newtonian case (see the dashed plot 
in the same figure) this ratio is in both relativistic cases 
for spinning and non-spinning black holes much larger in 
the inner accretion disk. The difference between the rel- 
ativistic and non-relativistic cases increases with the de- 
creasing height of the primary source of radiation. How- 
ever, we must recall that the incident radiation shown in 
the figures will still, after it is reprocessed by the disk, 
undergo some changes on its way from the disk to the 
observer caused by the relativistic effects. The emission 
from the vicinity of the black hole horizon will eventually 
be diminished by the redshift. As a result, for the low 

® Note, that we us e the power-law index a for the specific inten- 
sity, here, whereas in lDovciak et aLl l|2004bl ') power-law index F for 
the photon number density flux was used. 



lamp-post heights, the difference between the relativistic 
and Newtonian ratio of the received direct emission to 
the received refiected radiation will not be so large. 

4. REFLECTED EMISSION 

In the equatorial plane of the system we assume a neu- 
tral cold geometrically thin and optically thick Keple- 
rian accretion disk. Thus, neutral Fe fiuorescent Ka and 
K/3 lines are present in the local spectra together with 
the iron edge. Intensity of the reflected radiation in the 
local frame co-moving with the accretion disk was com- 
puted by the Monte Carlo multi-scattering code NOAR 
(jPumont et al.l[2000( ). Spectral features caused by scat- 
tering are hence automatically included in the model 
spectra, namely, the Compton hump is present; it occurs 
at energies typically ~ 15-50 keV (see example spectra 
in the left panel of Figure S]) . 

The normalization of the local flux as well as its shape 
(mainly in the Compton hump region) depend on the 
local scattering geometry, i.e., the incident and the emis- 
sion angles, as well as the relative azimuthal angle of the 
incident and emitted light rays. We stress the depen- 
dence on the incident and relative azimuthal angles of 
the lamp-post set-up. This feature is additional to the 
pure directionality dependence on the emission angle in 
the case of an extended corona when the disk is usually 
assumed to be illuminated isotropically. 

Th e single scattering approximation (jChandrasekhail 
|196C1| ) is used for the local polarization o f the reflected 
continuum component (for details see iDovciak et al.l 
'2004bD. The line flux and the primary radiation are sup- 
posed to be unpolarized. The polarization properties of 
the reflected radiation depend highly on the scattering 
geometry, as is illustrated in Figure |31 

The geometry of scattering is determined by the posi- 
tion where reprocessing happens on the disk. In this set- 
up the incident and emission angles are decided by the 
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Figure 3. Illumination of the disk. The radiation incident on the accretion disk is affected by the relativistic effects. The ratio of the 
incident radiation to the observed primary emission is shown for three heights of the primary source, h = 3, 15, and 100 GM/c^ (from 
left to right). The dependence on the radius is shown for the non- rotating Schwarzschild black hole (a=0) with the solid curve, for the 
maximally rotating Kerr black hole (a=l) with the dotted curve, and for the Newtonian case with the dashed curve. 




100 



0.6 - 



0.4 



0.2 







H3=0.0 ~ 




/ \ \ H-0-3 


H,=0.7 ~- 




/ \ \ npO.7 


He=0.9 - 









30 60 



90 

A* 



120 150 180 




60 90 120 150 180 

A* 



Figure 4. Reprocessed radiation. The properties of the local reprocessed radiation depend on the position on the disk. This is due to 
their dependence on the local scattering geometry. Left: The local intensity varies by more than one order of magnitude for different values 
of incident and emission angles. The shape of the Compton hump changes as well. Examples of spectra for three different pairs of values 
of the incident and emission cosines, /i; and fic, are shown. Here we assumed the relative azimuthal angle between incident and emitted 
light rays to be small (A<I> = 2.5°). Middle: The local polarization degree can have any value between and 1. We show the dependence 
on the relative azimuthal angle A<J> for three pairs of incident and emission cosines. Right: The local polarization angle obtains any value 
between —90° and 90°. The dependence on the relative azimuthal angle A<I> for three pairs of incident and emission cosines is shown. A 
polarization angle of 0° represents the direction perpendicular to the disk. 



incoming and outgoing geodesies. Throughout the disk 
all three angles, 5i,5c, and (Figure [T]), may attain al- 
most all possible values. Thus also the local polarization 
degree Pdisk and angle Xdisk nray achieve the whole range 
of values, to 1 and —90° to 90°, respectively, as is illus- 
trated in Figure ID In our computations, a polarization 
angle of 0° always represents the direction perpendicular 
to the disk and it is rotated counter-clockwise for positive 
values when looking towards the approaching photon. 

5. RADIATION FAR FROM SOURCE 

To get the resultant polarization far from the source, 
i.e. at radial infinity, one has to integrate the emission 
over the disk surface and thus one integrates the local 
polarization properties. From the previous section it is 
obvious that the result of such integration is not eas- 
ily estimated because of the complicated dependence on 
the local geometry of scattering. Moreover, relativistic 
effects enhance radiation from some parts of the disk. 
They also rotate the polarization angle, and thus the 
overall observed polarization properties will be given by 
an interplay of the local polarization properties on the 



disk and relativistic change acting on photons on their 
way to the observer. 

From our previous studies of po larization of therma. 1 
radiation of the accretion disk (see lDovciak et a"I]|2008D 
one expects relativistic depolarization of the local po- 
larization. This depolarization should be larger for a 
spinning black hole when the inner part of the accretion 
disk extends closer to the horizon, and thus is larger in 
size when compared to the Schwarzschild non-rotating 
case where there is a large central hole in the disk. It is 
in this inner part where the relativistic rotation of the 
polarization angle is the most variable. Yet, in the lamp- 
post geometry we may get different results (i.e. higher 
polarization for higher spin) when the primary radiation 
is taken into account as well, as we will see later on. 

The overall polarization at infinity is given by the 
Stokes parameters that have to be integrated over the 
accretion disk surface, 

-^o -^rcfl ^ -^prim 

= / d5G(/rk +/^k)+Vim, (4) 
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Figure 5. Reflected component at inflnity. The spectrum of the reflected component as the observer at inflnity would measure it for 
three different inchnation angles, 9o = 30° , 60° , and 80° (from left to right) and two heights of the primary source, h = 3 GM/c^ (solid) 
and h = 15 GM/c^ (dashed). The top panels correspond to the Schwarzschild case (a=0), the bottom panels correspond to the extremely 
rotating Kerr case (a=l). The grey lines depict the spectrum when the Fe Ka and K/3 lines are omitted. 



0.4 

0.3 

0.2 

0.1 


0.4 

0.3 

0.2 

0.1 







30° 








\ 

\ 


j 


1 


10 




\ 





e„=60° 



e„=80° 



0.2 



,2 0.1 



0.2 



,s 0.1 























X 





10 















A 















10 
E [keV] 



10 
E [keV] 




Figure 6. Polarization of the reflected component. The polarization of the reflected component as the observer at infinity would measure it 
for three different inclination angles, Qo = 30°, 60°, and 80° (from left to right) and two heights of the primary source, h = 3GM/c^ (solid) 
and h = 15 GM/c^ (dashed). The energy dependence of the polarization degree is shown. The top panels correspond to the Schwarzschild 
case (a=0), the bottom panels correspond to the extremely rotating Kerr case (a=l). The grey lines depict the polarization degree if the 
unpolarized Fe Ko and K/9 lines are not considered. 
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Here, the relativistic effects are expressed in terras of tfie 
transfer function G and the relativistic rotation of the 
polarization angle V'- For their definition and be havior 
see our previous work in [Dovciak (2004) ; Dovci ak et alJ 
(|2004al ). We have divided the intensity measured by the 
observer into the part coming directly from the primary 
source, /prim, and the part reflected by the accretion disk, 
/refl- The latter one can still be divided into the polar- 
ized continuum part, /^efl'' ^"^d the unpolarized intensity 



coming from fluorescent neutral iron lines (Ka and K/3), 
ll]^- We denote their counterparts in the local frames 
in which the radiation is emitted, i.e. before applying 
relativistic changes, as /lamp, Idisk, and /^"jj, respec- 
tively. The observed linear polarization degree Pq and 
polarization angle Xo are computed from the Stokes pa- 
rameters in the usual way, i.e. Pq = vQo + f^/-^o and 
tan2xo = Uo/Qo- 

6. RESULTS 

6.1. Polarization of the reprocessed radiation and iron 
line features 

First, we will look at the polarization properties of the 
relativistically broadened iron line. To this purpose we 




Figure 7. Reflected intensity at infinity. Left: The reflected component of tfie intensity at infinity integrated in the energy ranges 2-6, 
6-10, 10-20, and 20-50 keV (curves with lower to higher normalization) for three different observer inclinations, 8o = 30°, 60° and 80° 
(top to bottom panels). The dependence on the height of the primary source is shown by the solid curve for the non- rotating Schwarzschild 
black hole (a=0) and by the dotted curve for the maximally rotating Kerr black hole (a=l). Middle: The same as on the left but for the 
ratio of the reflected to the primary components of the intensity at infinity. Right: The same as on the left but for the dependence on the 
intensity of the primary source as an observer would measure it at infinity.In this graph, the height of the primary source increases with 
the primary intensity from to 100 GM/c^. 



show only the reflected component of the intensity /i-oa 
and polarization degree Picfi (Figures[5]and[6|). The value 
of the polarization angle does neither depend on the pri- 
mary radiation nor on the intensity of the iron lines, thus 
we will consider it only in the next section when we will 
discuss polarization properties with the primary emission 
taken into account as well. 

In Figure [5] one can see typical relativistic reflection 
spectra with a broad iron line below 10 keV and Compton 
hump that peaks around 20 keV. Here we show several 
examples from the lamp-post geometry varying different 
model parameters - the height of the primary source, 
the inclination of the observer, and the spin of the black 
hole. We show also the reflection continuum so that the 
intensity of the line is emphasized. The line is broader 
for lower h of the lamp-post because more photons hit 
the disk near its inner edge (see Figure [3]) where the line 
is shifted to the lower energy. 

Due to the fact that the Ka and K/3 lines are unpo- 
larized in our model, the polarization degree Piofi should 
be reduced at the energies to which the line is shifted. 
However, in some cases this needs not to be true and 
there even can be a rise in polarization degree, see the 
graphs for heights h = 15 GM/c^ in Kerr case in Fig- 
ure [6l The reason for this behavior is the presence of the 
iron edge, where the polarized continuum intensity sud- 
denly drops. As this edge is shifted to different energies 
there is a change in the contribution of the polarized flux 
at the energies near it. The overall effect can be either 
enhancing the polarization or decreasing it (see the grey 
graphs representing only the continuum emission). 

The polarization of the reflected component can be 
quite high. In our example figures it is the highest at 
lower energies, E < A keV, for lower heights of the pri- 
mary source, h — 3GM/c^, low inclination, 9o = 30°, 



and for the Schwarzschild black hole when it is slightly 
above 35%. 

Next we show the refiected emission (Figure [71) as a 
function of the lamp-post height, for the different en- 
ergy bands (2-6, 6-10, 10-20, and 20-50 keV), and for 
the three values of the inclination angle (30°, 60°, and 
80°). The reflected radiation first increases with the 
height, mainly due to the increase of the disk solid an- 
gle as seen by the source (i.e. the light-bending effect 
decreases). Thus for small heights the incident radiation 
on the disk increases with the height. The decrease at 
large heights (above h > 30 GM/c^) is primarily due to 
the finite size of the disk in our computations (we used 
Tout = 1000 GM/c^). However, even if the disk had ex- 
tended into infinity, there would still have been slight 
decrease in the refiected intensity with the height, as op- 
posed to the Newtonian case when it should stay con- 
stant. This is due to the fact that the reflected radiation 
from the inner disk is more intense in the relativistic case 
and as the lamp-post height increases the solid angle of 
the inner region, where relativity effects are important, 
decreases. 

For very large heights the reflected component ap- 
proaches the Newtonian case. From Figure [7] (left and 
middle panels) one can see that the intensity for differ- 
ent spins is almost the same already for heights around 
30 GM/c^, which is also due to the fact that at these 
heights the innermost accretion disk already makes quite 
a small contribution to the overall reflected radiation. 
From the figure it is also apparent that the contribution 
of the inner disk to the total intensity for the extreme 
Kerr black hole is the largest for high inclination angles. 
In Figure [7] we also show the dependence of the ratio of 
the refiected to the primary intensity on the lamp-post 
height and reflected versus primary intensity. Note, that 




Figure 8. Properties of tlie radiation at infinity. Intensity, polarized intensity, polarization degree, and polarization angle (top to bottom) 
at infinity for the total radiation, i.e. both primary and reflected components are included. The observer inclination is 60 = 30°, 60°, 
and 80° (from left to right) and the height of the primary source is h = ZGM/c^ (solid) and h = 15 GM/c^ (dashed). The black graphs 
correspond to the Schwarzschild case (a=0), the grey graphs correspond to the extremely rotating Kerr case (a=l). 



for low heights the contribution of the reflected compo- 
nent may be even larger than that of the primary one (in 
higher energy bands). All th ese effects are quali t ativel y 
similar to those described in iMiniutti fc FabianI ()2004D . 
who adopted a somewhat different geometry for the pri- 
mary source (i.e. a ring centered on the black hole axis). 

6.2. Overall polarization far from source 

If we consider the primary radiation as well, the de- 
crease of the polarization degree in the iron line region 
is much less visible, because the primary flux at these 
energies is still large. On the other hand, at high en- 
ergies, where the contribution from direct radiation is 
relatively lower and where the Compton scattering in- 
creases the contribution of the reflected radiation, the 
polarization degree increases. We illustrate this behav- 
ior in Figure [5] for two values of the lamp-post heights 
(ft, = 3 and 15 GM/c^). We also show the energy depen- 
dence of the polarization angle in the same figure. As 
previously mentioned for the polarization degree of the 
reflected component, one can see quite a large effect of 
the iron edge. 

The dependence of the polarization degree and angle 
at infinity on the lamp-post height, observer inclination, 
and total intensity at infinity are shown in Figures |9]- 



[T^ In all these figures both the reflected and primary 
emission are considered and the results are computed in 
several energy bands. 

To understand the relativistic results one has to take 
into account several aspects different from the Newtonian 
case: 

1. Local polarization. The emission angle, on which 
the local polarization properties depend, is not 
equal to the observer's inclination. Its value de- 
pends on the position on the disk from which the 
photon is emitted. The emission angle is approach- 
ing the observer inclination only for large radii. 
Therefore if the local polarization depended only on 
the emission angle, it would be the same all over the 
disk in the Newtonian case whereas it would change 
considerably in the relativistic case. The local in- 
cident angle and relative azimuthal angle between 
the incident and emitted light rays are different as 
well. 

2. Transfer function. The transfer function in the 
eqs. ([I])-® acts as a weight that amplifies the emis- 
sion in some areas on the disk while in other parts 
it can very efficiently suppress the outcoming ra- 




Figure 9. Polarization degree at infinity. Ttie values for the energy ranges 2 — 6, 6 — 10, 10 — 20, and 20 — 50keV (curves with lower 
to higher normalization) are shown. The solid curves are for the non- rotating Schwarzschild black hole (a=0), the dotted curves are for 
the maximally rotating Kerr black hole (a=l). Left: The dependence on the height of the primary source for three different observer 
inclinations, Bo = 30°, 60°, and 80° (top to bottom panels). Middle: The dependence on the observer inclination for three different heights 
of the primary source, h = 100, 15, and ZGMjf? (top to bottom panels). Right: The dependence on the total intensity at infinity for three 
different observer inclinations, 6o = 30°, 60°, and 80° (top to bottom panels). In this graph, the height of the primary source increases 
with the total intensity from to lOOGAf/c^. 



diation. Due to the relativistic effects this impact 
is non-axisymmetric. Thus whereas in the New- 
tonian case the overaU polarization angle may be 
parallel or perpendicular due to the symmetry of 
the system, this does not apply to the relativistic 
case where the observed polarization angle may, in 
principle, attain any value. 

3. Energy dependence. While in the Newtonian case 
the energy dependence of the local polarization 
does not play any important role it is not so in 
the relativistic case. Due to energy shift of the pho- 
tons the local intensities in the eqs. (|l])-(|n]) work as 
weight similarly as the transfer function mentioned 
above. Thus the features such as the unpolarized 
iron line, the iron edge and the Compton hump 
may change the contribution of different parts of 
the disk to the polarization in a particular energy 
band, especially if those features are shifted in and 
out of the energy band of interest. 

4. The dependence of the relativistic rotation of the 
polarization angle on observer inclination. To un- 
derstand this effect one has to realize the role of 
the critical point, where the photons are emitted 
perpendicularly to the disk. The existence of this 
point is due to special relativistic aberration. For 
small inclinations even very small orbital speed of 
the matter moving towards the observer is enough 
to ensure that light rays emitted perpendicularly 
to the disk reach the observer. For very large in- 
clinations the orbital speed must be much higher. 
Therefore the critical point is far away from the 
center for small observer inclinations whereas it 
moves close towards the center for very high ones. 
In the latter case the contribution of the general 



relativistic aberration and light bending becomes 
important and as a result the location of the criti- 
cal point is behind the black hole. For small incli- 
nations the critical point is at the azimuth where 
the disk is moving towards the observer. For the 
importance of the special relativistic effects (aber- 
ration, Doppler effect and thus special relativistic 
rotation of the polarization angle) and the gen- 
eral relativistic effects (significant for small radii: 
light bending, gravitational shift, and general rela- 
ti vistic rotation o f polar ization angle) see Figure 3 
in 'D ovciak et al.l ()2008[ ) , Appendix D in iDovciakI 
(20041), or Figure 1.3 in iDovciakI (|2010l ). O ne can 
see there that for radii below the critical radius 
the relativistic rotation of the polarization angle, 
-0, spans the whole possible range (i.e. from —180° 
to -f 180°) whereas above the critical radius the ro- 
tation is restricted to some interval that gets quite 
narrow for large radii. 

5. Spin of the black hole. The disk around the black 
hole is supposed to terminate at the marginally sta- 
ble orbit which means there is a cavity in the disk 
around the black hole; the lower the black hole spin, 
the larger the cavity. Because the largest influence 
of the relativistic effects are visible mainly in close 
vicinity to the black hole, the relativistic effects are 
much more pronounced for highly spinning black 
holes. 

As a consequence, different parts of the disk will emit 
light with different local polarization. The direction of 
this polarization will be further changed as the photon 
travels from local frame co-moving with the disk to the 
stationary observer and this change will be again depen- 
dent on the place of the emission. The contribution to 
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Figure 10. Dependence of the polarization angle at infinity on the height of the primary source. A polarization angle of 0° represents 
the direction perpendicular to the disk. The values for the energy ranges 2 — 6, 6 — 10, 10 — 20, and 20 — 50 keV for three different observer 
inclinations, Qo = 30°, 60°, and 80° (top to bottom panels) are shown. Left: The non-rotating Schwarzschild black hole case (a=0). Right: 
The maximally rotating Kerr black hole case (a=l). 



the overall polarization front each part of the disk will 
be different depending on the weight given by the trans- 
fer function and the local intensity (both dependent on 
the position on the disk). The latter one is proportional 
to the incident intensity from the primary source thus 
the height of the lamp determines at which radii the disk 
shines the most. The incident angle under which photons 
strike the surface of the disk also depends on the height 
of the lamp, so the height modifies the local polarization 
degree and angle as well. 

Thus while in the Newtonian case it is quite straight- 
forward to see how the polarization by reflection from 
the disk behaves it is much more complicated in the rel- 
ativistic case. Therefore, in the following, we will restrict 
ourselves to the description of the most interesting results 
and we will discuss only the main qualitative features. 

We begin with a bit surprising and contra-intuitive re- 
sult — the polarization degree is highest when the radia- 
tion from the inner accretion disk has a large contribution 
to the overall flux. This is opposite to the de-polarization 
one gets in th e inner disk for the Comptonized thermal 
radiation, (see lDovciak et aP 2008: S chnittman fc KrolikI 
[20091 The cause of the thermal de-polarization is the 
relativistic rotation of the polarization angle. The local 
polarization is in this case always either parallel or per- 
pendicular to the disk. In the reflection model the local 



polarization angle depends on the local scattering geom- 
etry. When integrating the contributions from the disk, 
the local polarization angle adds up with the relativistic 
rotation in such a way that the overall polarization is 
much less de-polarized. 

This has two consequences. Firstly, the highest polar- 
ization degrees are achieved for low lamp-post heights 
when the inner region of the disk is more illuminated 
than its more distant parts. Secondly, the radiation from 
the extreme Kerr black-hole accretion disk is more po- 
larized than that in the Schwarzschild case where there 
is a hole in the disk below marginally stable orbit. Note, 
however, that there are some exceptions to the latter be- 
havior, see the results for some intervals of the lamp-post 
heights and for the observer inclination 30° and 80° . This 
can be explained by the fact that the polarization in the 
Schwarzschild case is determined only by the disk above 
marginally stable orbit, r > r^ns — 6GM/c^, while in 
the Kerr case the contribution from the disk below and 
above this radius have different polarizations, effectively 
reducing the overall polarization. 

Note, that in Kerr case the dependence of the polar- 
ization degree on the height. Figure |9l roughly follows 
the dependence of the ratio of the reflected and primary 
radiation, /lofl/^prim, see middle panel in Figure [T] This 
is not true for the Schwarzschild black hole where for low 
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Figure 11. Dependence of the polarization angle at infinity on the observer inclination. A polarization angle of 0° represents the direction 
perpendicular to the disk. The values for the energy ranges 2-6, 6-10, 10—20, and 20-50 keV for three different heights of the primary 
source, h = 100, 15, and 3GM/c^ (top to bottom panels) are shown. Left: The non-rotating Schwarzschild black hole case (a=0). Right: 
The maximally rotating Kerr black hole case (a=l). 



lamp-post heights the polarization is quite large in spite 
of the fact that the reflected intensity is much smaller 
(due to the missing innermost part of the disk). This, 
again, is just another consequence of the fact that the 
contribution of the inner disk above the marginally stable 
orbit to the overall observed polarization is significant. 

We should also mention that similarly to the intensity 
also the polarization degree and angle do not change with 
the spin for large heights {h > 30 GM/c^) of the primary 
source (see left panel in Figure [5] and both panels in Fig- 
ure [lOl). 

The polarization degree has a maximum for low ob- 
server inclinations, e.g. in the highest studied energy 
band, 20-50 keV, for the lamp-post height h = Z GM/c^ 
it peaks in the extreme Kerr case at inclination 25° where 
it reaches almost 15% polarization (see the middle panel 
in Figure 19]). 

The polarization degree for zero inclination angle has 
to be zero because of the symmetry. As soon as the sym- 
metry is broken by non-zero inclination the total polar- 
ization increases, determined by the region below critical 
radius which for low inclinations is far from the center. 
As soon as the emission above this radius starts to con- 
tribute significantly the overall polarization starts to de- 
crease with the inclination angle (the larger the inclina- 
tion the lower the critical radius) . This turn-around is at 



lower inclinations for larger heights because higher lamp- 
post illuminates better farther radii (i.e. radii above the 
critical one). There is another turn-around at the incli- 
nation when the critical radius moves close to the cen- 
ter. The contribution to the polarization from below and 
above critical radius cancel each other and the depen- 
dence of the polarization degree on the observer inclina- 
tion reaches its minimum. For even higher inclinations 
the polarization again increases, it is determined mainly 
by the emission from far above the critical radius with 
lower and lower contamination from the regions around 
and below this radius. For very high inclinations the 
reflection is small and thus the overall polarization de- 
creases again. 

The polarization angle at infinity is quite sensitive to 
the details described in the points mentioned above and 
its dependence on height and inclination is rather com- 
plex. Note mainly the different behavior in different en- 
ergy bands (Figures fTUHTT]) . While the polarization de- 
gree behaves in different energy bands very similarly, just 
being scaled up with higher energy, the polarization an- 
gle dependences may differ substantially. As discussed 
above for the polarization degree, the polarization angle 
for low inclinations is also determined mainly by the re- 
gion below critical radius (top panels in Figure [TU]) while 
for high inclinations by the region above this radius (mid- 
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Figure 12. Dependence of the polarization angle at infinity on the total observed intensity. A polarization angle of 0° represents 
the direction perpendicular to the disk. The values for the energy ranges 2-6, 6—10, 10—20, and 20-50 keV for three different observer 
inclinations, 9o = 30°, 60°, and 80° (top to bottom panels) are shown. The height of the primary source increases with the total intensity 
from to lOOGAf/c^. Left: The non-rotating Schwarzschild black hole case (a=0). Right: The maximally rotating Kerr black hole case 
(a=l). 



die and bottom panels in the same figure) . The change in 
behavior is very well visible by rapid change of polariza- 
tion angle with inclination in Figure 1111 The transition 
happens at the inclination when the polarization above 
critical radius starts to dominate for higher inclinations. 
The transition depends on the lamp-post height and en- 
ergy band and it can be either gradual or quite abrupt. 
Notice that there is no transition (or a very mild one) in 
the bottom panel in Figure [TT] for the extreme Kerr black 
hole and for a very low height of the primary source. This 
is due to the fact that the inner region below the critical 
point still has large impact on the polarization at infin- 
ity even for very high inclinations. In this case also the 
dip in the dependence of the polarization degree on in- 
clination is not so deep, see the dotted graphs of bottom 
middle panel in Figure [S] 

For very high inclinations the polarization angle at in- 
finity is determined by outer regions of the disk. Here 
relativistic effects are small enough so that the polariza- 
tion at infinity is given by the local values in the disk. 
Because the incident and emission angles are very large in 
this area and for this inclination, the polarization angle 
will be 0° (or 180°). The values in Figure [TT] are smaller, 
mainly for high heights because the disk in our computa- 
tions docs not extend to infinity but up to 1000GA//c^ 
and therefore the incident angle is smaller than 90° and 



the local polarization angle decreases. 

It is good to understand the behavior of the polariza- 
tion for different lamp-post heights and observer inclina- 
tions, however, the height cannot be directly measured 
and to test the dependence on inclination we would need 
polarimetric observations of many AGNs. Therefore in 
the next sections we will use the dependence of the po- 
larization on the directly measurable total intensity that 
can be observed at infinity. To this purpose we show 
the theoretical curves for different energy bands for the 
polarization degree in the right panel in Figure [3] and 
for the polarization angle in Figure [121 In these figures 
the height of the primary source increases with the total 
intensity. 

7. OBSERVATIONAL PERSPECTIVES 

7.1. Monte Carlo simulations of the detector response 

The complex dependence of the polarization on the 
geometry of the system gives important insight in the 
emission process at work in accreting black holes. Many 
missions with the capability of measuring polarization 
are expected or are discussed for launch in the next few 
years and therefore it is important to investigate if the 
sensitivity is sufficient to perform polarimetry also of rel- 
atively faint sources like AGN. For this purpose, we de- 
veloped a dedicated Monte Carlo software which is de- 
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scribed below. 

Instruments able to measure polarization exploit the 
modulation which appears in the case of polarized radia- 
tion in the azimuthal response. Experimentally, this re- 
sponse, or modulation curve, is the histogram of the ini- 
tial d irection of photoelectron for photoelectric polarime- 
ters (jBellazzini fc Spandrel . |2010() or that of the s catter - 
ing direction for Compton devices (jMcConnelll . 120100 . 
If the impinging photons are polarized, the photoelec- 
trons (or the scattered photons in the case of Compton 
instruments) are not isotropically distributed, but their 
azimuthal distribution on the plane orthogonal to the 
direction of incidence shows a cosine square modulation. 
The peak (or the minimum) corresponds to the direction 
of polarization, while the amplitude of the modulation is 
linearly proportional to the degree of polarization. The 
maximum modulation achieved for 100% polarized radi- 
ation is called modulation factor of the instrument /i and 
it is in general a function of the energy. 

In the following we will assume instruments based 
on a phot oelectric polarimeter, the Gas Pixel Detec- 
tor fGPD: iBeiTazzini fc Mulenl I20T0I ). The systematic 
effects of GPP have been reported to be below 0.3% 
(|Bellazzini &: Spandrel l2010l ). and therefore they can be 
neglected if the polarization degree is above ~1%. In 
this hypothesis, the sensitivity of a polarimeter is only 
a matter of how many photons are collected. The num- 
ber of photoelectrons emitted in each azimuthal angu- 
lar bin is Poisson distributed and its fluctuations ulti- 
mately limit the capability to measure the amplitude of 
the modulation curve a nd, consequently, the polarization 
(IWeisskopf et al.llMol) . 

The basic idea behind our Monte Carlo software is to 
derive the degree and the angle of polarization by gener- 
ating test modulation curves, that is trial histograms rep- 
resenting the azimuthal distribution of photoelectrons. 
The entries in each histogram, one for each spectral in- 
terval of interest, is the number of collected photons in 
that energy band. The number of counts Ni in each 
angular bin i is distributed as a cosine square function, 
whose phase and amplitude are the angle of polariza- 
tion and the product of the modulation factor and the 
degree of polarization averaged in the energy band, re- 
spectively. A Poisson noise -s/ZVi is added in each angu- 
lar bin. The modulation curve is fitted with a function 
M{(j)) = A + B cos{(f> — (j)o), where A and B are con- 
stants. The fit returns the "measured" angle of polar- 
ization = (po and its error a^j, while the polarization 
degree Vj (alongside with its error) is derived from the 
usual formula 

max ~f~ -^^min 

fi B + 2A' 

The fit is repeated J times, e.g. 1000, adding a different 
Poisson noise contribution to the histogram for each trial. 
Eventually the expected ("measured") degree and angle 
of polarization (and their errors) are derived by averaging 
the values obtained by the j-th trial, 

7^=^= j^T'j, av^^v, CT* = CT$,. (8) 

j 

This method was used by iDovciak et al.l (|2008( ) to de- 
rive the sensitivity of next X-ray polarimetry missions 



to the rotation of the plane of polarization in the emis- 
sion from accreting Galactic black holes. In that case 
polarization was to be measured in quite narrow energy 
bands, about 1 keV, and therefore it was possible to de- 
fine the average quantities required to generate the test 
modulation curves, e.g. the modulation factor, the an- 
gle and the degree of polarization. Instead here we need 
to study the sensitivity on a wider energy interval and 
therefore we present a different approach, initially pro- 
posed by Costa (private communication), which has also 
the advantage of taking into account the finite energy 
resolution of the instrument. 

The simulation starts from the knowledge of the mod- 
ulation factor, the spectrum of the source, the expected 
degree and the expected angle of polarization. They are 
calculated in K narrow spectral bins, typically 0.1 keV 
which is well below the energy resolution of the Gas Pixel 
Detector (about 1 keV at 6 keV). The values of the mod- 
ulation factor, of the spectrum, of the degree and of the 
angle of polarization in the fc-th spectral bin Ek will be 
named fj,k, Fk, Pk, ^Pk, respectively. For each /c-th bin in 
energy, we build a histogram with a cosine square mod- 
ulation proportional to Pk fik and a phase tpk ■ The his- 
tograms have / angular bins between and 2tt, e.g. 72 
bins corresponding to intervals of 5 degrees, and the to- 
tal number of elements is proportional to Fk- The value 
of the fc-th element in the i-th angular bin Hki is basi- 
cally the fraction of the flux in the spectral bin Ek which 
causes a photoelectron to be reconstructed in the i-th 
angular direction (neglecting statistical fluctuations) . So 
far this is not very different from the approach formerly 
presented, with the only difference that we define his- 
tograms in very narrow energy bands. 

The finite spectral resolution of the instrument makes 
the measurement of the energy of the photon to be 
smeared around an average value. This causes a pho- 
ton at energy E to be reconstructed as a photon at dif- 
ferent energy E, and the difference between E and E 
depends on the energy resolution. The uncertainty on 
the energy measurement needs to be convolved with the 
response of the instrument, in particular with the depen- 
dence of the modulation factor on energy. As a matter 
of fact, completely polarized photons at energy E cause 
a modulation of amplitude n{E) which can be quite dif- 
ferent from fJ-{E). To model properly the amplitude of 
the modulation, one should take into account that the 
instrument responds to polarization differently with the 
energy and that the energy is measured with a finite res- 
olution. In our Monte Carlo simulations we include this 
degeneracy. We assume, and this is appropriate in most 
cases, that the effect of the finite energy resolution is 
limited to smear the energy of the event and that it does 
not act on the direction of emission of the photoelectron, 
namely on the polarization. Therefore we apply the blur- 
ring due to the finite spectral resolution to the spectra 
resolved in angular bins. Practically, this is done by tak- 
ing the K values Hf^j for any angular bin i = i and for 
k G [0,K — 1], where we recall that k is the index running 
over the spectral bins. These / spectra are introduced 
in XSPEC as additional models and the response matrix 
of the instrument is used to apply a Gaussian smearing 
to them with the XSPEC tool fake it. The results are 
/ pha files which contain the number of counts Nki col- 



14 



Dovciak, Muleri, Goosmann, Karas & Matt 



0.10 



0.00 





Figure 13. Correlation of the degree (left) and of the angle of polarization (right) with the flux for MCG-6-30-15. Each point corresponds 
to a different height of the source, specified within the figure (in units of GMjr?). It is assumed that each state of the source is observed 
for 500 ks with the Medium Energy Polarimeter on-board NHXM in the 8-30 keV energy interval. The dashed and the dotted lines are 
the expected dependency for spin 1 and spin 0, respectively. 



lected in the angular bin i for each spectral bin in the 
detector energy space. 

The arrays N^n for any fixed value of fc' and any value 
of i represent a sort of modulation curves defined in 
very narrow energy bands E^i of 0.1 keV that we name 
"pseudo modulation curve" (because the Poisson noise 
has not been included yet). Any practical measurement 
will be performed in much larger energy intervals, of the 
order of a few keV. Therefore we build the pseudo modu- 
lation curves in any energy band of interest lS.Ei by sum- 
ming those histograms N^i which are referred to spectral 
bins contained in the interval Ai?;. This results in L 
pseudo modulation curves, as many as the number of in- 
tervals, which we name M;, where we drop the subscript 
i (we implicitly assume the angular dependence being 
indicated by i). The values of Mi take into account the 
effect of the finite energy resolution on the amplitude of 
the modulation and that the degree and the angle of po- 
larization can change with energy, that is these quantities 
are assumed constant only on small intervals of 0.1 keV. 

After that the pseudo modulation curves are generated, 
the Monte Carlo scheme proceeds as discussed above. 
In each lS.Ei energy band, the software produces J trial 
modulation curves, which are calculated from the same 
pseudo modulation curves Mi but adding for each trial 
a different Poisson noise. Each trial modulation curve 
is fitted with the function M.{<^) = A + B cos(0 — 0o), 
and eventually the "measured" degree and angle of polar- 
ization are derived by averaging the values obtained by 
each trial; see equation ([8]). The value of the modulation 
factor used in each energy band to derive the polariza- 
tion degree from the amplitude of the modulation (see 
equation [7]) is obtained by weighting the dependency of 
fi with energy on the spectrum of the source. 

7.2. Results 

The Monte Carlo scheme was used to investigate if po- 
larimeters on-board next missions could detect distinc- 
tive signatures of the lamp-post model. A key character- 
istic of such a model is that the height of the illuminating 
source must change to explain the temporal variation of 
the source. This results also in a change of the observed 



fiux and polarization, as reported in right panel in Fig- 
ure [9] and 1121 Therefore our primary objective is to dis- 
cuss if next missions will be able to find an evidence of 
the positive correlation between the degree of polariza- 
tion and the fiux. 

7.2.1. Predictions for AGN case: MCG-6-30-15 

First we model the case of Seyfert 1 galaxy MCG-6- 
30-15. From the simulations presented above we know 
how the flux and the polarization correlate, but we still 
need to associate the primary source height to a total flux 
measured in a certain state. Since we can not derive this 
normalization from the available inform ation, we assume 
that the lowest flux state measured by iLee et ahl (|20000 
corresponds to h = 3GM/c^. This choice, though arbi- 
trary, is supported by conclusions of other authors, who 
suggest that the illuminating source is quite close to the 
central black hole. We assume that the observation in- 
clination is 30° and that the black hole in MCG-6-30-15 
is maximally rotating. 

An interesting result expected on the basis of the lamp- 
post model is that the polarization usually increases with 
energy (see Figure [8]). Therefore in the following we 
will take as an example the New Hard X-ray Mission 
(NHXM hereafter), which is dedicated to perform imag- 
ing, spectroscopy, and polarimetry in a wide energy inter- 
val ext ending from soft to hard X-rays (jTagliaferri et al. l. 
I2010al) . NHXM will allow for polarimetry above 10 keV 
where the reflected radiation start to dominate the emis- 
sion. In the focus of one of the four multilayer telescopes 
there will be alternatively placed two photoelectric po- 
larimeters based on the Gas Pixel Detector, which are 
sensitive in the intervals 2-10 and 6-35 keV, and pos- 
sibly a third Compton polarime ter which wou ld extend 
the response up to 80 keV (Soffit ta~et al.ll20l"0l ). Here we 
focus on the Medium Energy Polarimeter (MEP), the 
polarimeter that will be sensitive between 6 and 35 keV. 

We show the correlations of the degree and of the angle 
of polarization with the fiux in Figure [T31 Both the flux 
and the polarization are integrated in the energy range 8- 
30 keV to avoid the iron line which is expected to be un- 
polarized. Each point represents an observation of 500 ks 
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Figure 14. The same as Figure [Tsl but for the Galactic black hole system XTE J1650-500. It is assumed that each state of the source is 
observed for 50 ks with the MEP on-board NHXM in the 8-30 keV energy interval. 



and corresponds to a different state of the source, that 
is a different height as indicated in the figure. Ahhough 
the angle of polarization is almost constant with the flux, 
a certain evolution of the degree of polarization should 
be detectable. Moreover, we report for comparison the 
evolution of the polarization expected in the case of a 
Schwarzschild black hole to show that such evolution is 
somehow dependent on the spin of the central black hole. 

7.2.2. Predictions for Galactic black hole systems: XTE 
J 1650-500 

Features which can be regarded as due to reflection, 
such as relativistic iron lines, are ob served i n a n umber 
of Galactic black hole systems; see iMilleii ()2007D for a 
review. There have been attempts to apply the light 
bending scenario also to these sources in analogy with the 
AGN case and therefore we discuss also the possibility to 
observe the correlation between the polarization and the 
flux in the case of Galactic objects. 

We take as an example the transient black hole sys- 
tem XTE J1 650-500 and assoc iate the maximum flux 
measured by iRossi et aD (|2005l ) to the a source height 
h = %GM/c^, changing the source height between this 
value and h — "iGMjc?. We assume that the source 
spin is 1 and the observer inclination is 30°. We sup- 
pose to observe each state of XTE J1650-500 for 50 ks 
with the MEP on-board NHXM. The resuh is reported 
in Figure [TH In this case a significant detection could be 
reached with quite short observations. 

8. DISCUSSION AND CONCLUSIONS 

In this paper we have discussed the observational prop- 
erties of the intensity of the emerging primary radiation 
component, the radiation component incident on the ac- 
cretion disk, and the reflected radiation in the lamp-post 
geometry (the light bending model). We have extended 
our previous work on the subject and showed new results 
for polarization computations, including the unpolar- 
ized iron line complex in 2-10 keV energy range and the 
Compton hump at higher energy band (above 15 keV). 
We used our theoretical computations to model possi- 
ble future observations by next-generation X-ray satellite 
missions equipped with an X-ray polarimeter on board. 



For low inclinations we find that the polarization of 
the reflected radiation is diminished in the broad iron 
line energy region. For high inclinations, however, the 
situation is more complicated due to the shift of the iron 
edge; this effect enhances the observed polarization and 
it may balance the expected decrease. In principle the 
effect of lower polarization of the relativistic broad line 
region could be used to discriminate between the rela- 
tivistically broadened iron line and the partial covering 
scenario, which are both proposed to explain the excess 
of flux in the spectra of black-hole accret ion disks be- 
tween 2-10 keV (see e.g. lMiller eralll2009[) . Note, how- 
ever, that we did not try to estimate the behavior of the 
polarization in the partial covering scenario in this paper. 

The polarization degree in a lamp-post geometry is 
higher at the highest studied energy band, 20-50 keV. 
This is a natural result coming from the fact that the pri- 
mary source has a power-law spectrum with a negative 
index. Thus the best results for polarimetry of reflec- 
tion spectra should be achieved in the Compton hump 
energy region. We also get higher polarization degree for 
large values of spin of the black hole, small height of the 
primary source, and low inclination of the observer. 

As an example of the polarimetric sensitivity of next 
missions, we could assume the performance of t he Grav- 
ity and E xtreme Magnetism SMEX (GEMS; iJahodal 
2010: Blac k et al.l . [201Q) . which will be launched by 2014. 
However, the energy range of the polarimeters on-board 
GEMS will be limited between 2 and 10 keV, with the 
peak of the sensitivity being around 3 — 4 keV. Thus it 
will observe photons in a region where the iron (unpo- 
larized) fluorescent line is prominent and the direct ra- 
diation from the corona largely dominates the reflected 
component. The reflected radiation is the only one po- 
larized in our model. Consequently the expected de- 
gree of polarization is quite low, between 1 and 2% 
depending on the height of the illuminating source. 
Thus the signal could be more easily affected by even 
a small amount of polarization of some other origin, 
e.g. by Comptonize d thermal emission in case of X-ray 
Black hole binaries (iDovciak et al.l [20081: iLi et al.l [20091: 
ISchnittman fc KrolikL I2009ir Therefore we preferred to 
simulate observations of the polarization in higher energy 
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band. 

In our model there were several assumptions and ap- 
proximations that could affect the resulting polarization. 
First of all the primary radiation was unpolarized. It is 
clear that the effect this will have on our results will be 
most prominent in low energy bands, for higher lamp- 
post heights, low inclinations or small black-hole spins 
where primary emission is dominant (see Figure [T]). Thus 
the polarization at low energies could be enhanced but 
determined by the primary source while the polarization 
at higher energies could be reduced. 

In our computations we used the single scattering ap- 
proximation for local polarization properties. The accu- 
racy of this approximation depends on the energy. One 
can have an idea of the effect confronting Figure 1 (lower 
panel) and Figure 3 in Matt (1993). It follows from these 
figures that for energy below lOkcV results should not 
change much, while they could be different by about 25% 
in the 10 - 30 keV and by almost 40% in the 30 - 50 keV 
energy range. 

We assumed a smooth equatorial disk. In case of a 
"rough" disk, with random ripples of the length much 
smaller than the length at which relativistic effects 
change considerably, the local polarization properties 
should be averaged over incident and emission angles. 
Because we average over these values anyway when in- 
tegrating over the disk but including the relativistic ef- 
fects it is hard to predict in a unique way how this will 
change the overall polarization at infinity. In case the 
disk is smooth but not strictly equatorial (warped disk 
but still close to the equatorial plane), the local geom- 
etry of scattering in the disk will be changed and thus 
the polarization at infinity will be changed as well, both 
dependent on the geometry of warping. In case of largely 
warped disks or non-aligned disks we expect different re- 
sults mainly because the disk would be more illuminated 
at places closer to and facing the primary source, and 
because of the different local geometry of scattering and 
different values of the relativistic e ffects. 

Our work is complementary to iSchnittman fc KrolikI 
(|20Tnh who assume a different geometry and structure of 
the corona. A direct comparison of their paper with our 
results is not easy because of different approximations 
adopted in these two works and the notorious sensitivity 
of polarization to details of the model. Nevertheless, the 
general set-up is similar in both scenarios and the level 
of expected polarization is comparable. 

The behavior of the polarization degree and polariza- 
tion angle in the light bending model is quite complex. 
Because the result depends on the interplay of several 
parameters, the polarization properties may be degener- 
ate with respect to different parameter values (e.g. for 
larger heights we cannot distinguish a difference between 
a rotating and a non-rotating black hole). Nevertheless, 
when combined with spectral and timing observations, 
polarimetry is yet another important channel that can 
help us to uncover the physical parameters of the black- 
hole accretion disk systems, such as the black hole spin, 
the system inclination, and the height of the illuminating 
source. 
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